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The problem being addressed 


Casualty calculations for any given orbit ASSUME that an object enters 
with equal probability over all possible locations along an orbit arc. 
(Equally distributed in Argument of Latitude (“ArgLat” or O) 


Simulations and qualitative physics show that the equatorial bulge of 
the Earth pushes a “wall of air” into the path of the object at the 
equator, equal to several times the density scale height, and potentially 
up to a 50x local density increase at the equator relative to the density 
at the poles. 


Periodic geometric and gravitational variations are both distorting the 
orbit 


THEREFORE: 
— there should be periodicity! 
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Why it is Cyclic 


(and why it precedes the nodal crossing) 
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Why it is not simple: 
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Final Decay Parameters vs. Latitude Sy 


¢ The equatorial de =e 
bulge affects 
drag. Note the 
concave altitude 
and density plots 
vs. latitude. 


¢ However, there 
are interesting 
(and countering) ne | rn i ee : 
J3 convex effects 
going on in the 
radius plot 
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Prior Work 


Matney and Bacon had independently hypothesized an effect due to the oblate 
Earth 


—  Matney, M. J., “Statistical Issues for Uncontrolled Reentry Hazards,” 3rd IAASS Conference, Rome, Italy, 2008. 
— Bacon, J.B. “The Clustering of Natural Decays in Equatorial Zones” Internal ISS Whitepaper 8/17/2015 
Matney published an early study of real polar decays, indicating a clustering of 


decays near the equator. 


—  Matney, M. J., “Statistical Issues for Uncontrolled Reentry Hazards: Empirical Tests of the Predicted Footprint 
for Uncontrolled Satellite Reentry Hazards”, 4th IAASS Conference, Versailles-Paris, France, 2011. 


Bacon and Matney collaborated on an IAASS workshop paper in May 2016 
confirming that the effect was significant, demonstrating a sinusoidal 


compression and rarefaction of decays, and suggesting forward work. 


— Bacon, J.B. & Matney, M.J. (2016). Statistical Issues for Calculating Reentry Hazards, IAASS 8th Workshop on 
Space Safety, May 21 2016. 


Bacon at Darmstadt 2017 provided additional insight into the mechanisms and 


structure of the decay compression 


— Bacon, J.B. (2017). Minimum dV for Targeted Spacecraft Disposal Proc. 7th European Conference on Space 
Debris, Darmstadt, Germany, 18-21 April 2017 


Claire Fremeaux and her group at CNES suggested in Darmstadt 2017 that a J3 
gravitational effect was also at work, evident in her decay studies of highly 


eccentric orbits. 


— Cid Borobia, E., Fremeaux, C., and Goester J.F.: “Fast Re-Entry Deorbitation with Acceptable Risk 
Level” Proc. 7th European Conference on Space Debris, Darmstadt, Germany, 18-21 April 2017 
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So...The New Approach 


We go looking for periodicity (cyclic effects) in decay profiles 
— Both symmetric and asymmetric around the equator 
— Convert analysis space to Argument of Latitude, not Latitude 


¢ Normalizes to common reference curve 

¢ Directly models the concentration effect in ways meaningful to Ec 

¢ Latitude can be reconstructed: ArcSin[{Sin(Inclination)}*{Sin(Argument 
of Latitude)}] 


To be complete, we need to study a LOT of profiles 


To characterize periodicity, we formulate the clustering effect as 
a Fourier series, and look for the major periodic terms. 
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First we simulate a lot of entries: 


We ran the General Mission Analysis Tool (GMAT) on 
over 30K decay cases to 90 km: 

— 18 inclinations from 5-90 in 5-degree steps 

— 4 Ballistic numbers (50, 100, 150, 200 kg/m2) 

— 4Right Ascensions of Ascending Node at the Vernal Equinox 

— 101 incremental BN changes to smoothly dither the 4-day decay 

¢ AND, there are about 18,000 propagations per run, for > 0.5E9 propagations! 

Several dozen follow-up cases were also run for 
comparison: 
101 dithered runs each for: 

— Entry on Summer solstice at 100 BN, Inclination =60, 24 RAANs 

— BN=600 kg/m2 all inclinations, Vernal Equinox RAAN = 90 

— BN= 300, 400, 500 at 90 inclination Vernal Equinox RAAN = 90 
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100 entries around Arg(Lat) for Every Case Na ASA 


..for every Inclination, BN, and RAAN 


First, find how many points (N points) lie along 360 degrees for 
Argument of Latitude 
— Typically not an integer: e.g., 57.632 


The average expected gap in Argument of Latitude to the next point is 
Average = 360/N 
— Inthe example, Average = [360/(57.632)] = 6.247 degrees 


Measure the actual gap in Argument of Latitude to the next point 


Ratio the actual gap to the average gap 
— Inregions of clustered decays, the number is less than one 
— Inregions of rarified decays, the number is greater than one. 


PLOT the compression ratio vs Argument of Latitude 
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Comression Value (Dimensionless) 


Sample Compression File 
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Fourier Transform the Compression Curve 


Divide Argument of Latitude into 256 even increments 
Sort the actual compression curve points by Argument of Latitude 


LOOK UP the closest Compression Curve value to each of the 256 
Argument of Latitude sample bins 


— NOTE: we only have 100 points, and thus a chunky/noisy binned 
representation, BUT we only need the 1° 4 terms (long averages) 


Fourier transform the 2 
bin-sampled curve te sa! sr ST 


— (XCEL application) 


Collect amplitude and phase of first four Fourier terms 
— Constant (Zeroth term) is forced equal to ONE. 
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Compare Sizes and Phases of the Fourier 
Terms 


¢ Even-numbered terms are symmetric North-South 


(32, J4) 


— The Cos(20) term is dominant (30-40% compression factor) 


¢ It varies slightly with RAAN and in a sine-like way with inclination 

¢ This is due principally to the “wall of air”, plus perhaps a J2 Gravitation effect 

¢ We don’t yet understand the “S-curve” growth of amplitude with inclination* 
The Cos(40) term is 1-3% compression factor 

¢ This generally rises with inclination, from near zero 


e Odd-numbered terms show North-South biasing 


The Cos(O) term is highly erratic and variable with inclination and 
RAAN, generally < 10% and >1% compression 

¢ We believe that this is solar heating dominating a small J3 effect 
The Cos(30) term is 1-3% compression and is very repeatable over 
all conditions. 
The Cos(5O0) term is negligible, and not expected from J4 Gravity model 
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RAAN Family of Compression Curves: 


Solstice/ Inc=60 / BN=100 / RAAN_Variable Compression Curves 
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@ RAAN=240 @ RAAN= 255 —@—RAAN= 270 @ RAAN=285 @ RAAN=300 @ RAAN=315 @ RAAN= 330 —@—RAAN= 345 
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Average 10 Fourier Term Compression Effect (Dimensionless) 


The J3 Term Falls Out When Averaged! 


Averaged Effect of First (10) Fourier Term Over All RAAN's at 60 Inclination 
0.03 


° When the first Fourier term in 10 is 
averaged at each value of © over all 
RAANS, a perfect sine wave 
emerges, representing the residue 
believed to be the J3 gravitational 
contribution. This slightly biases 

natural decays to occur in the 

northern (vastly more populated) 
hemisphere. (2.7% bias) 


0.02 


0.01 


-0.01 


-0.02 


-0.03 


0 30 60 50 120 150 180 210 240 270 300 
Argument of Latitude © (Degrees) 


14/32 


330 


360 


JBB 


~< | 
AG 
(o> 
& 1 
y| 


d and Fourth Fourier Amplitudes 


Sulat; 
UU 


Secon 


0.35 


I 
I 
! 
I 
I 
| 
\ 
\ 


- 


oO 
4 
Zz 
< 
oO 


RAAN 90 


RAAN 0 


\ 
\ 
i} 
1 
t 
1 
I 
! 

U 


U 


~ 
~ 
~ 


=-s==s_ 
—— - 


ee 


Ue 
ee 


oe 
EPSARSANRN een 


~s 


bod 


0.3 


~ 


NN 


0.05 BNSO “.,BN100 BN150 B 


7 


N200 


i. 


N200 
\ 
} 


N200 
i] 
i 


aseiaAy T JEUIWON punole apnyjdwy sauno4 


7 


687 


LvZ 
Tre 


677 


ST 
6eT 


Sequence Number (Arbitrary) 


2Amp #@4Amp 


JBB 


15/32 


, Third, & Fifth Fourier Amplitudes 
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Compression Curve First 


National Aeronautics and Space Administration 


0.35 
0.3 
0.25 
0.2 
0.15 
0.1 


aSe1aAy T JEUIWON punoue apnyjdwy sauno4 


— —_— —_—— 9A 
es 


= Jha 


S67 
687 
c 
L£L7 
TLZ 
S9T 
657 
eGC 
Lvé 
Tee 
SEC 
677 
ECC 
L£TZ 
TTc 
SOC 
66T 
e6T 
£8T 
T8T 
SLT 
69T 
e9T 
ZST 
TST 
SvT 
6eT 
eet 
£eT 


STT 
60T 
e0T 


Sequence Number (Arbitrary) 


@iAmp #@3Amp #5Amp 


JBB 


16/32 


National Aeronautics and Space Administration 
Compression Curve Components at 90 
Inclination 


Combined-Even and -Odd Fourier Terms High Inclination 
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Compression Curve Components at 10 
Inclination 


Combined-Even and -Odd Fourier Terms Low Inclination 
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Amplitude and Phase of Fourier Terms vs. BN 


First Five Fourier Amplitudes With Changing BN First Five Fourier Phases With Changing BN 


700 
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First Five Fourier Amplitudes With Changing BN 


20 term dominates amplitude 
All other terms a <5% effect 


Only the even terms show monotonic phase 
behavior as BN varies 
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To Apply Decay Compression Factor to E;: 


e We now have an analytic curve of how concentrated or rarified the 
entries are around any given orbit, vs. Argument of Latitude 
¢ We incrementally step N (>>1) steps in Argument of Latitude to 
make a 360 arc. 
¢ From our analytic representation, we have the concentration 
factor (C(O)) for that portion of the arc 
Remember: C(O) is < 1 for locally concentrated decays! 
¢ We calculate the latitude at that particular step /; 
Latitude = ArcSin[{Sin(Inclination)}*{Sin(Argument of Latitude,)}] 
¢ We look up the average population in that latitude band 
— Note: we get to the same latitude twice, with different compression 
factors at each passage!) 
¢ We DIVIDE the average population per km? by the compression 
factor C(O) 


e (Note! A binomial correction is needed!) 


¢ We sum the result in each step, and divide by N 
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So: What we do: 


We have a dumb model: 
— Find the net risk increase from a summary risk chart of the already-calculated 
cases. 
¢ Because of the huge computation effort, needs the assumption of reasonably 
uniformly-distributed net population growth by latitude band (not true..). 
Ora slightly better model 
— (parametric fit published in the paper...) 


Or an even better model: 


— Linearly interpolate phases and amplitudes between those of the bracketing 2 
inclinations and BNs (a big spreadsheet...), and apply this best-estimate curve 
to exact latitude band population forecast. 


Ora best model: 


— Just run an automated GMAT series at all RAANs for the BN and inclination 
desired, and then apply the averaged Fourier curve to the best population 
model 


« ESPECIALLY if we know the decay season or date. 
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Crude Model (in Paper) 


We set out to makea 
universal, simple model 
(orange curves) of the 
averaged Fourier terms 
(blue curves) 


4 cosine amplitudes, 4 
phases 


Some terms were easy to 
fit simply, some weren’t 


A lot of interpolation may 
miss fine details 
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Conclusions: 


¢ The wall-of-air is a real effect that grows in importance as 
inclination and BN increase 


¢ Diurnal (10) effects are real and large, but average out 
over the seasons and orbit precessions. 
— These may be important if one knows exact entry date 


¢ The higher-order symmetric (40) and asymmetric (30) 
terms are purely gravitational, and about 1/6 as important 
as the atmospheric terms 


— Although small, because the factors are highly repeatable over all 
BNs and RAANSs they can be added as real biases to the 
compression curves 
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And most importantly: 


¢ We recommend that analysts use a Fourier model with 
only the 20, 30, and 40 derived terms of a specific GMAT 


entry simulation* that will represent the average entry 
concentration. 


¢ The compression curve for the specific case should then 
be applied to the population bands as per slide 21. 


— (*Simulation conditions include specific inclination and BN of the 
spacecraft to 80 km + estimated footprint length downrange as 
Argument of Latitude, at the Vernal equinox, 270 RAAN) 
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BACKUP 
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True confessions 


These curves don’t take us to the ground. We need to get 
there analytically before calculating Ec. 


The method described does not account for lateral 
spread of the footprint 


We have not proven yet that the 10 term averages to zero 
— In work, see hot-off-the-presses chart 


We don’t have physics explanations (yet) for some of the 
interesting knees/inflections in the Fourier constant 
trends 
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Potential Issues (so far...) 


There’s no such thing as a circular orbit or perfect initial conditions 


1. BUT: Minimum dV entry study and current runs show that the entry is driven mostly 
by nonlinear growth of deceleration in the last 3 orbits, and is thus very insensitive 
to initial conditions three days earlier. Chose to make orbits approximately equal 
altitude at ascending and descending node. 


Decay to 90 km speeds up the runs by nearly a factor of four vs. 
propagation to 80. 


¢ Bacon study in April 2017 showed little change in remaining downrange gap from 90 to 
80km over all BNs and inclinations, but this adds some fuzz 


¢ For E, we need to calculate to the GROUND, and must add in a correction factor 


Running intentionally at (near) the Vernal equinox eliminates a 
suspected diurnal bulge which can put the sub-solar point up to 23.5 
degrees off the equator as a major density perturbation 

¢ Weare just starting to see if this subsolar location effect averages out over a year. 

¢ See Hot-off-the-presses chart 

The spread over 1, 2, or often 3 orbits of dither shows MINOR changes 
in the compression curves. These are smoothed and averaged, but 
not characterized. 
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Potential Issues (so far...) Part Il 


5. The Fourier binning of the data adds some scatter 
vr am The premise is that this does not matter in the first few terms, whose period is much 
~ a greater than the raw data gap 


6. Differing precession rates and decay times mean that there was some 
small variation in the RAAN (and solar vector) location of the entry: it 
only STARTED at 90-degree increment RAANs 


7. The compression curve algorithm only generates gaps to the next 
point forward. To be mathematically correct, one should average with 
the gap backwards too. 


8. Sampling using VLOOKUP in EXCEL gets the closest value BELOW 
the one you are looking up. There is no telling whether the next higher 
one is 1000x closer, without complicating every equation. 
¢ This affects the Fourier binning process. 
¢ Both this and the prior issue affect the phase of the compression curve slightly 
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Potential Issues (so far...) Part Ill 


9. The latitude equation ArcSin(Sin(Inclination)*Sin(ArgLat)) 
does not account for the elliptical geoid effects. 


10. The future population in any given latitude band is stilla 
bit of a guess, despite Chris Ostrom’s recent excellent 
work! 


11. The compression curves are inverted to some people’s 
liking 
° It made intuitive mathematical sense to normalize each gap value to some 
local constant (the average gap width) rather than the inverse. This getsa 
curve in normalized gap, whose inverse at any local point is the 
“compression” 


° Unfortunately, the AVG(1/X) is not 1/AVG(X), and the resulting curves need a 
small tweak to renormalize the ultimate PDF to equal 1 over all latitudes. 
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J3 Effect 


e From Fremeaux et al 2017: 


— Indeed, due to the Earth potential term J3, re-entry orbits with 
arguments of perigee in the North hemisphere are more prone to a 
faster descent than the orbits with South perigee arguments, as shown 
in Figur e 3. Re-entry duration 
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— Figure 3. Re-entry duration depending on perigee argument (w) and solar activity. 
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